Kekule-distortion-induced Exciton instability in graphene 
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Effects of a Kekule distortion on exciton instability in single-layer graphene are discussed. In the 
framework of quantum electrodynamics the mass of the electron generated dynamically is worked 
out using a Schwinger-Dyson equation. For homogeneous lattice distortion it is shown that the gen- 
erated mass is independent of the amplitude of the lattice distortion at the one-loop approximation. 
Formation of excitons induced by the homogeneous Kekule distortion could appear only through 
direct dependence of the lattice distortion. 
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I. INTRODUCTION 

Recent experimental accessibility of graphenei^ has 
drawn much interest on this material which present a 
wide variety of interesting properties^iSiZiS. Most of the 
properties of the graphene arises from the peculiar energy 
spectrum near the so-called Dirac nodal points 9 -*^. 

Recent work on exciton instability in graphene 
monolayers is based on the Dirac Hamiltonian 
descriptio n 11 ' 12 ! 13 . The exciton gap is derived and 
solved through a self-consistent equation similar to 
the one appearing in the chiral symmetry breaking 
phenomenon^ 4 . It was shown that an exciton can 
be formed under a strong long-ranged particle-hole 
interactio n 12 : 15 . Exciton can also be formed in a single- 
layer graphene through the mechanism of magnetic 
catalysis of dynamical mass generation as pointed out 
ir*l£. This work showed that the magnetic catalysis can 
induce exciton condensation even for weak particle-hole 
coupling 15 . These results are obtained in the framework 
of quantum electrodynamics QED deduced from the 
linear energy spectrum of the graphene monolayer. 

Exciton instability in graphene bilayer systems have 
been studied in the case of a short-ranged Coulomb 
interaction and a finite voltage difference between the 
layers-1 Self-consistent exciton gap equations are de- 
rived in the framework of Hartree-Fock approximation 
and it is shown that a critical strength of the Coulomb 
interaction exists for the formation of excitons. The crit- 
ical strength depends on the amount of voltage difference 
between the layers and on the inter-layer hopping param- 
eter. The voltage difference drives a gap in the energy 
spectrum of the graphene bilaye r 18 : 19 and combined to a 
strong Coulomb interaction leads to an exciton instabil- 
ity. 

Similarly a gap can open in the energy spectrum of 
graphene monolayers by means of a lattice distortion 
such as the so called Kekule distortio n 20 : 21 ' 22 . A weakly 
screened Coulomb interaction and a gap in the linear en- 
ergy spectrum of the graphene monolayer are favorable 
elements for the formation of excitons. 

We focus our attention here on the consequences of a 
Kekule distortion on dynamical mass generation for elec- 
trons in the graphene monolayer. Indeed quasiparticle in 



2D systems can acquire mass through dynamical interac- 
tion between electrons and holes leading to an increase 
of the energy gap 14 : 23 ' 24 . We address the question : can 
a Kekule distortion affect the dynamical mass generation 
mechanism in such a way that the resulting energy gap 
favors the formation of excitons ? We show that in the 
specific case of an homogeneous Kekule distortion and at 
the one-loop approximation the dynamical mass genera- 
tion is unaffected by the lattice distortion. We show that 
the gap in the energy spectrum is a sum of two indepen- 
dent contributions : one induced by the homogeneous 
Kekule distortion and one induced by dynamical mass 
generation. Consequently exciton instability can only be 
formed through direct dependence of the amplitude of the 
lattice distortion and the Kekule-independent dynamical 
mass of the electrons. There is no amplification effects of 
the homogeneous Kekule distortion by the mechanism of 
dynamical mass generation. 

The outline of the paper is the following. In section |TT] 
we recall the construction of the quantum electrodynamic 
action of the graphene monolayer in presence of a Kekule 
distortion. In section lllTl we present the calculation of the 
dynamical mass term using the Schwinger-Dyson equa- 
tion of the electron. Section HVl summarizes and discusses 
the present work. 



II. QED 3 ACTION OF THE GRAPHENE 
MONOLAYER WITH A KEKULE DISTORTION 

A graphene monolayer is a honeycomb array of atoms 
of carbon as depicted in Fig. [TJ In a monolayer the elec- 
trons can hop between nearest-neighbour carbon atoms 
through 7r-orbitals with energy t. The unit cell of the 
graphene monolayer is composed of two type of carbon 
atoms and we denote them by A and B. The whole lat- 
tice is devided in to two sublattices, the A^-type and 
the A^-type. The Hamiltonian describing the graphene 
monolayer reads 



H a = -t 



E 

vGAa a— 1 



^2 alb 



r+r Q 



(i) 



2 




FIG. 1: Representation of the graphene monolayer with the 
nearest-neighbour vectors r a connecting the two sublattices 
A a and A^. 



where a and b stand for the 7r-electron creation and ani- 
hilation fermionic operator on the atoms of type A and B 
respectively. The vector r a connects the two sublattices 
A a and A# and is given by r a = ae a with a the lat- 
tice parameter and e\ = (1,0), e 2 = (— l/2,v / 3/2) and 
e 3 = (— 1/2, — V3/2) are unit vectors. In the following 
the lattice parameter a will be set equal to one and plays 
the role of the unit of length. 

The diagonalisation of the Hamiltonian dXJ) leads to 

the kinetic energy Ek = t\ J2 a =i e th ' ec " |, where e Q 's are 
the nearest-neighbour vectors of the graphene monolayer. 
The kinetic energy vanishes at the two independent nodal 
points which are chosen as K + = ^0, t^tj) an d 

K- = —K + in the Brillouin zone. At low energy the 
bare Hamiltonian H$ of the 7r-electron can be rewritten 
by considering the energetic contributions arround the 
nodal points and reads^ 



-J2 v Fvl(r){-2d z )v b {r) + h.c. 



where we made use of the notations in the complex 
plan z = x + iy and d z = i (d x — id y ). The velocity 



vp = ^ta/h where we set h 



1 (the kinetic en- 



ergy becomes e p = hvp\p\). The amplitude u a (b) and 
v a (b) are smooth functional operators and are connected 
to the creation and annihilation operator a and b through 
the Fourier components u a (jp) = a P +K + ,Ub(p) — b p+ K + 
and v a (p) = a p+ K_,Vb(p) = b p+ K_ where p is the wave 
vector. The Hamiltonian Ho describes the hopping of the 
7r-electrons in an undistorted lattice. 

A lattice distortion can be modeled by a variation of 



the hopping parameter 8t r ^ a which depends on the posi- 
tion in the lattice r and on the direction e Q . The corre- 
sponding Hamiltonian for the graphene monolayer reads 



H K = - X! X! 5t r,aa\.b r+ba + h.c. 



(2) 



In the following we consider the Kekule distortion for 
which the bonds between the carbon atoms are arranged 
similar to the benzene molecule 20 . The lattice is then pic- 
tured by an alternation of long and short bonds between 
the carbon atoms of types A and B^. The Kekule dis- 
tortion is modeled by the hopping parameter 8t rtQ , and 
reads St r , a = § (A(r)e lK + r «e lG - r + A(r)e ,Jf - r °e llGT ), 
where A(r) stands for the amplitude of the Kekule dis- 
tortion and are the nodal points in the Brillouin 
zone. The vector G = K + — K- connects the two inde- 
pendent Dirac cones located at the nodal points K + ^_y 
Focusing on the low-energy contributions of the Kekule 
distortion around the nodal points the Hamiltonian Hk 
is given by 



H K = - J2 [ A ( r ) u K r ) v b ( r ) + A(r)»t (r)u 6 (r)] + h.c. 

r£A A 

The Hamiltonian describing the 7r-electrons in the 
graphene single-layer with a Kekule distortion is then the 
sum of the bare and Kekule Hamiltonians H = Hq + Hk ■ 
Using the fermionic spinor — [ui,(r)u a (r)v a (r)vi,(r)]^ 
the Hamiltonian H can be rewritten in the quantum elec- 
trodynamic framework and reads 



H = / d 2 rifj(r) 



VFlkdk + A(r) 



A(r) 



The aim of the present 



where the gamma matrices are defined by 70 = T3 ® T3, 
7i = Vpn (8 t 3 and 72 = vpT 2 (8 t 3 , T{ 01 2 } are the Pauli 

matrices. The element A(r) is a 4 x 4 matrices related 
to the Kekule distortion amplitude A(r) by the relation 
A(r)r 3 
-A(r)r 3 

work being to characterize the behaviour of the dynam- 
ical mass generation for graphene monolayer in presence 
of an homogeneous Kekule distortion we reduce our study 
to the case A(r) = Aq. Finally for a system at temper- 
ature T — 1/(3 the action of the graphene monolayer is 
given by 



Spi = 



dr Id ftp(r, r) 



7^ + A r 



V(r,r). (3) 



The presence of the electromagnetic field surrounding 
the graphene monolayer is also to be considered. Indeed 
the graphene monolayer is embedded in an electromag- 
netic field which spans the whole 3D space. However 
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the electrons are confined in the 2D space delimited by 
the carbon atoms. Hence the density of charge and cur- 
rent verify p(x,y,z) = 8(z)p 2 D{x,y) and J(x,y,z) = 
S(z)J2D{x,y). The electromagnetic field vector a and 
scalar <p potentials are related to the Green functions 

of the laplacian 2v — V 2 rather than V 2 for the three- 
dimensionnal electromagnetic fiel d 23 ' 26 . To describe the 
electromagnetic field embedding the graphenc monolayer 
it is convenient to use the following Euclidean QED 3 ac- 
tion 



III. DYNAMICAL MASS GENERATION 



A. The photon propagator at finite temperature 



Integrating over the fermion fields ip leads to a pure 
gauge Lagrangian C a = ^a^A^a^ where A M „ is the 
dressed photon propagator from which we shall extract 
an effective interaction potential between two fermion 
and derive the dynamical mass of the fermions. 

The finite temperature photon propagator in Eu- 
clidean space verifies the Dyson equatio n 29 ' 30 



d A x- 



1 



2V z a 2 



ffJ,L> f 

J J {ivi 



(4) 



where = d^a^ — d^a^ is the electromagnetic field ten- 
sor derived from the three-dimensionnal vector potential 
a M . The Fourier transform of action ((4} leads to the bare 

photon propagator A^°J = (8^ — Q^Qu) /(2\q\) in Eu- 
clidean space. For finite temperature the imaginary-time 
component qo of the (2+l)-dimensionnal wave-vector 
Q = (qo,qi,q2) is given by the bosonic Matsubara fre- 
quency qo — 2im/ j3 (n is an integer and takes its values 
in the range ] — oo, oo[). 

Finally the full quantum electrodynamic action de- 
scribing the graphene monolayer with an homogeneous 
Kekule distortion and embedded in a 3D electromagnetic 
field reads S = S e .m. + S e i 



S = 



d 3 G 



-1 



2\/^F 



+ / d 3 xip 



7 M {dp. - iga^) + A 



tf>. 



(5) 



It is known that 2D systems described by a quantum 
electrodynamic action like ([5]) experience a dynamical 
mass generation for the fermionic field ip in interaction 
with an U{\) gauge field. Appelquist et ah 14 ' 24 showed 
that at zero temperature the originally massless fermions 
can acquire a dynamically generated mass when the num- 
ber N of fermion flavors is lower than the critical value 
N c = 32/ V 2 . Later Maris 2 - confirmed the existence of a 
critical value N c ~ 3.3 below which the dynamical mass 
can be generated. Since we consider only spin-1/2 sys- 
tems, N — 2 and hence N < N c . At finite temper- 
ature Dorey and Mavromatos22. and Lee2£ showed that 
the dynamically generated mass vanishes at a temper- 
ature T larger than the critical one T c . More recents 
works have been performed on dynamical mass genera- 
tion in graphene monolaye r 11 ' 12 . However the question 
now arise about the effects of a lattice distortion like 
the Kekule distortion. How does the mass generated dy- 
namically in presence of a Kekule distortion behave? We 
concentrate on the effects of an homogeneous Kekule dis- 
tortion and for small temperature T — > 0. 



A" 1 



n 



/if 7 



where the bare photon propagator A^°J is derived from 
the action of the bare electromagnetic field ^ and the 
polarisation function is obtained from the integration 
over the fermionic field tp in the action (|5|). 

For the computation of the dynamically generated 
mass it is enough to consider the static temporal 
dressed photon propagator component A no {q° = 0, q) 
for which the bare photon propagator component reads 
A Q ° \q° = 0,q) = l/2\q\. The detailed calculation of the 
static temporal component of the polarisation function 
rioo(<Z = 0, q) is given in appendix lAl and reads 



U°°(q = 0,q) = 




In (2 cosh (7r9 g )) 



(6) 



where (3 is the inverse temperature and 9 = 
(j^j ~ ^v^q 1 + l^ol 2 - The coupling parameter 

a is related to the number of fermion flavor N = 2 and 
to the electron charge g, a = Ag 2 N. 

For very small temperature T — > the polarization 



function 



becomes 



n u V = o,0 



8nv 2 F 
(v F qY 



|A | 
4|A | 2 



2v F \q\ 



arctan 



VFq 
2|An| 



The process of dynamical mass generation is dominated 
by mechanisms at large wave vectors q. Indeed for q ^> 
2 1 Aq\/vf the polarisation function is asymptoticaly equal 
to ts- 2 — I q I which confirms the fact that the dynamics 
between electrons and holes are dominated by large wave 
vectors 1 ^. 
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B. The electron self-energy 

We now derives the electron self-energy which is also 
the dynamical mass. The Sch winger-Dyson equation for 
the electron propagator at hnite temperature reads 



G^~\k) 



d 2 p 



(27T): 



:l ti G{p)A tlv {k -p)T v , 



(7) 



where p = (po = ujF,n>P), G is the dressed electron prop- 
agator, r„ the electron-photon vertex which will be ap- 
proximated here by its bare value gj v and A^ is the 
dressed photon propagator. The second term in ([7|) is 

the fermion self-energy £, (G- 1 = G(°) - E). Per- 
forming the trace over the 7 matrices and working out 
the sum over the fermionic Matsubara frequencies u>F,n 
in equation ([7]) leads to a self-consistent equation for the 
self-energy of the form 



£(fc) = g 1 



d 2 p 
{2-Kf 



A O0 (0,fc 



p) }^ t a nil 



where ep 



.(0)2 



2£p 



ef- + E(p) 2 and ef 2 = v F f + |A | 2 are 
respectively the energy spectrum of the graphene mono- 
layer with and without correction of the one-loop ap- 
proximation near the nodal points K + t_y In the limit 
q 3> 2|Ao|/u_f the dressed photon propagator is given by 
Aoo(0, q) = [(2 + a/(W7TVF))\q\]~ ■ The angular integra- 
tion in equation ([5]) is achieved using the approximation 
f(\p-k\) = 0(p-k)f(\p\)+9(k-p)f(\k\) where / is 
a function which depends on the absolute value of its 



j3ep 
2 



(8) 



arguments 
reduces to 



13.15 



Hence the self-consistent equation (|8]) 



E(jfe) = C 



dp- 



pE(p) 



dp 



k\/p 2 + {A(p)/v F ) 

S(p) 
p 2 + (A( P )/v F f 



(9) 



: Ag 2 / (32tto f + a), 
E(p) 2 and the ultraviolet 



where we introduced the constant C 
the function A(p) = ^\A \ 2 
cutoff A. The electron self-energy is assumed to be a 
small correction to the electron propagator. It follows 
that for k ^> \Ao\/vf the denominator of the second 
term in equation ^ is approximately equal to p. In the 
first term of equation the function A(p) plays the role 
of a cutoff for p — > and equation © can be rewritten 



dp^+j-dp^i 



S(p) 



)A(p=0)/v F 



. A derivation 



of this equation up to the second order with respect to 
the wave vector k of the electron self-energy leads to the 
differential equation 



fc 2 E (k) + 2/cE (k) + CE(fc) = 0, 



(10) 



for which the infrared and ultraviolet boundary condi- 
tions are given respectively by 



fc 2 E'(fc) 
(fcE'(fe) + E(fc) 



k=A/v f 



fc=A 



= 0, and 



= 0. 



(11) 
(12) 



The differential equation (fT0|) admits a solution which 
verifies also the boundary conditions (fTTjl and (f!2|) and 
readsHii^ 



£(fc) = 



A 3 / 2 /tan(<5) , f v F k 

; sm In — = — 

sin((5) v / ^i V 2 \ A 



(13) 



where the phase 8 is equal to arctan \J AC — 1. In or- 
der to verify the boundary condition (|12[) the ampli- 
tude A(p = 0) has to verify the following relation A = 
(v F A) exp [-45/V4C - lj . 

The solution of the electron self-energy has been de- 
rived for wave vectors larger than the cutoff \Aq\/vf and 
shows to be independent of the amplitude A(r) = A of 
the homogeneous Kekule distortion. The independence 
of E on the amplitude Ao implies that the dynamical 
mass generation mechanism cannot be controlled by an 
homogeneous Kekule distortion. 

It was shown elsewhere in the context of graphene bi- 
layer that a gap opening in the energy spectrum can lead 
to an exciton instability under the specific amplitude of a 
short-ranged Coulomb interaction 17 . Here a similar exci- 
tonic instability could appear since the graphene mono- 
layer spectrum sees a gap open directly induced by the 
Kekule distortion. Moreover the gap is enlarged by for- 
mation of a mass which is dynamically generated. An ex- 
citon instability could take placed 3 .. We address the ques- 
tion whether a Kekule distortion can amplify the mecha- 
nism of dynamical mass generation and consequently af- 
fect the exciton instability. The solution (jT2J) for the elec- 
tron self-energy E shows that the dynamically generated 
mass of the electron is independent of the homogeneous 
Kekule distortion. However the mass of the electron 
given by m(k) = \/|Ao| 2 + E(fc) 2 = (f (fc)#(fc)) shows 
that the exciton instability resulting from the electron- 
hole interaction is directly dependent on the homoge- 
neous Kekule distortion but not through the mechanism 
of dynamical mass generation. 
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IV. CONCLUSIONS 

The electron self-energy for graphenc monolayer with 
an homogeneous Kekule distortion has been derived in 
the framework of quantum electrodynamics. It has been 
shown that the dynamical mass generation of the elec- 
trons resulting from the electron-hole interaction in the 
graphene monolayer is independent of the amplitude of 
the homogeneous Kekule distortion for a one-loop ap- 
proximation. Such an independence of the dynamically 
generated mass provides an insight on the effects implied 
by lattice distortion on the energy spectrum of graphene 
monolayer. 

The gap in the energy spectrum of the graphene mono- 
layer is relate d to the mass o f the electron which is equal 
to m(k) = y/\A \ 2 + E(fc) 2 = (V(k)V(k)). The Kekule 
distortion control directly the gap by means of the ampli- 
tude |Ao| without any amplification through the dynam- 
ical mass £ (due to the independence mentioned previ- 
ously). Consequently exciton instability can be formed 
by direct relation to the homogeneous Kekule distortion 
and by the Kekule-independent mechanism of dynamical 
mass generation. 

So far we tried to shed light on the effects induced by 
lattice distortions such as an homogeneous Kekule distor- 
tion A(r) = A . The present study does not admit any 
conclusion on the independence of the dynamical mass 
m{k) in presence of an inhomogeneous Kekule distortion 
A(r) ^ const. It is expected that inhomogeneous lat- 
tice distortion affects significantly the energy spectrum of 
graphene monolayer. As a consequence exciton instabil- 
ity could possibly be amplified through lattice distortion 
via non-trivial dynamical mass generation. 



gamma matrices. Focusing only on the static temporal 
component of the polarisation n o(<z° = 0, q) and using 



the following relations Tr 



iA o7 °iAo7 c 



-4|A | 2 and 



Tr [Yk p .7°7 v (k v + q v h°] = ig { p ^k p (k v + q v ) where we 

defined the metric tensor g^ 1 ^ — diag(l, —v F , —v F ), the 
polarisation function can be reduced to 



<J = ± Ui v 



d 2 k / 1 



fc 2 + |A | 



I 



(fc + g) 2 + |Ao| 



x { 4gWk p (k n + q v )~4\ A \ 2 



(Al) 



The computation of the polarisation function can 
be simplified using the Feynmann identity — 

Jo dx (ax+(l-x)bf ' applying the change of variable k -> 

2tt 



' (ax + (l-x)b) 2 

k —xq (in other words w'f = '^S- (n + 1/2) and k' = k+xq) 



and performing the sum over the fermionic Matsubara 
frequencies ujj one gets 




dx{ Si -2 



(v F k 2 



\A Q \ 2 )S 2 
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APPENDIX A: DERIVATION OF THE 
POLARISATION FUNCTION 

The polarisation function rLj„ is obtained by integrat- 
ing the QED 3 action ([5]) over the fermionic field ip and 
reads 



We define by Si and S' 2 the sums over the Matsubara 
frequencies given by the following relations^ 



Si 



So 



E 



n— — oo 



L' 2 + v 2 F k' 2 + x(l - xjvlq 2 + |A | S 



4ttF 

OO 

E 



tanh(7rF) 



1 



„ ... [w' 2 + v 2 F k' 2 + x(l - xjvlq 2 + |A | 2 
f3 2 1 dSi 



IP"(g) 



9_ 





EE 



d 2 k 



Ti[G (k)Y-G (k + q)Y] 



where Go(fc) = -pqli^pr is the electron Green function 
and the trace operator Tr runs over the space of the polarisation function reads 



where u> = *£-(n + 1/2). The integration over the wave 
vector k can be performed through the change of variable 
Y = 4;\/ v Fk' 2 + X ( X ~ x ) v f^ + l A o| 2 and finally the 
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n°V = o,a) 



— dx lim / 

P Jo A ^°°./e c 



2tt 



ft + r 



<9Y 



dx lim 



A^oo J A 27T 



x(l — x)^ 



ay 



dX \2vvlP 



ln(2cosh(7re 9 )) 



/3|A |y tanh(7r6 g ) 



where g = 



(£) ^(l-^'f + IAol 2 
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